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Abstract. Many cosmological models of gamma-ray bursts (GRBs) assume that a 
single relativistic shell carries kinetic energy away from the source and later converts 
it into gamma rays, perhaps by interactions with the interstellar medium or by shocks 
within the shell. Although such models are able to reproduce general trends in GRB 
time histories, it is difficult to reproduce the high degree of variability often seen in 
GRBs. We investigate methods of achieving this variability using a simplified external 
shock model. Since our model emphasizes geometric and statistical considerations, 
£> ■ rather than the detailed physics of the shell, it is applicable to any theory that relies 

on relativistic shells. We find that the variability in GRBs gives strong clues to the 
efficiency with which the shell converts its kinetic energy into gamma rays. 

(N 

The "external shock" models of gamma-ray bursts (GRBs) assume that the 



gamma rays are emitted from a single shell of material traveling at highly rela- 
tivistic speeds (for example, [1] and [2]). Fenimore, Madras, and Nayakshin [3] 
have found that the envelope of emission for a single relativistic shell fits the over- 



all shape of a few GRBs, but it does not account for the wide variability found in 
GRBs. In this paper, we investigate methods of achieving the observed degree of 



variability in GRB time histories using randomly placed active regions on the shell. 



"3 : DESCRIPTION OF THE MODEL 



The "time" that is measured in GRB time histories is the time at which photons 
arrive at the detector (denoted by T), and is not the time of emission as would 
be measured in the rest frame of the detector (equivalent to the rest frame of the 
shell's expansion, denoted by t). We set t = to be the time at which the central 
explosion occurs, and T = to be the arrival time of a photon emitted at t = 0; 
therefore 

T = (1-/5 cos 6)t = A 1 - 1 t (1) 

where A is defined as A = 7(1 — j3[i), and \i = cos 9. By setting Equation 1 equal 
to a constant, one finds that the surface of constant arrival time is an ellipsoid 
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FIGURE 1. a. (left) The position of the spherical shell at four different times as it expands 
from t to i max . The two ellipses represent surfaces of constant arrival time. Photons emitted 
from any point on the inner ellipsoid will arrive at the detector simultaneously at time T, and 
photons emitted from any point on the outer ellipsoid will arrive at the detector at time T + dT . 
The cross-hatched region represents the volume that will be "seen" by the detector between these 
two times. Each small circle represents the volume swept out by an entity during its emission 
lifetime of AT p . Only the white entities will be seen between T and T + AT. b. (right) A 
three-dimensional view of the cross-hatched area shown at left. The volume of this annulus is 
given by T(T)dT. 



[4]. All of the photons emitted from such an ellipsoid will arrive at the detector 
simultaneously. 

Our model considers a thin shell (as required by the observations [3]) that ex- 
pands outward from a central point so that R = [3d,. The shell expands to a radius 
of R = (3ct , and then emits gamma rays until it reaches a radius of -R max = /3c£ max . 
The gamma ray emission must occur from small, independent patches on the shell, 
which we call "entities" [3]. Each entity begins as a small perturbation, perhaps 
caused by the shell's interactions with the interstellar medium or by instabilities 
within the shell. The entity grows at the sound speed c s ~ c until it reaches a 
maximum size of AR± = TcAT p (see Table 2 in [3]). Thus, each entity represents 
a causally connected region, and many entities can fit on the surface of the shell. 
We assume that each entity emits isotropically with a power-law photon spectrum 
in its own reference frame, $ (E ) = E ~ a photons (entity dt dE dVL ) _1 (primed 
quantities are measured in the rest frame of the entity). 

Figure 1 gives a pictorial representation of our model. It is important to empha- 
size that we are considering a thin shell. Thus, the four dashed lines in Figure la 
represent the shell at four different times. The area between t Q and t max represents 
the volume swept up by the shell over its emission lifetime. Likewise, the entities 
are small areas on the shell. The small spheres in Figure la represent the volumes 
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swept out by the entities as the shell expands. 

A more quantitative description can be made by calculating the volume enclosed 
in the cross-hatched region, which we denote as T(T)dT (Figure lb): 

/■min(R max ,R ell out ) [-^(T+dT) 

T(T)dT = 2tt / / r 2 dp dr (2) 

Jnim(R ,R cll ,in) •> »( T ) 

where -R c ii,out = l3c{T + dT)/(\ — (3), and -R e ii,m = /3cT/(l + (5). Equation 1 gives 
fi(T) = [ft- 1 -Tc/r], and p(T + dT) = [/T 1 - (T + AT)c/r}. Evaluating the 
integral yields the volume seen per dT: 

T(T) = if T < T 

nc (/3c) 

= (l-/3) 2 

TIC (/9c) 

= (l-/3) 2 

= j^ (tL, - r4(1 T ^ )4 ) if r 2 (i + /3) 2 T < T < T 2 (l + /?) 2 T max 

= ifT>r 2 (l + /3) 2 T max 

where T Q — t /(l — /3), and T max = £ max /(l — 0)- Note the rather surprising result 
that the volume "seen" by the detector is constant for T max < T < T 2 (l + [3) 2 T for 
the relativistically expanding shell. During these times, the number of entities in 
the cross-hatched region in Figure la is a constant with respect to T. In comparison, 
in the non-relativistic case, the volume seen, and therefore the number of entities 
seen, would increase as the area of the shell increased. 
The time history from such a shell is given by 

/•mm(R max ,_R cll . out ) r-^{T+dT) 

V(T)dT = 2tt / / pC'A^ a+1) r 2 dp dr (4) 

where V(T) is the expected time history in photons (dT dAd e t)~ l , C is the pho- 
ton flux as observed in the rest frame of the entity, and hr^ a+1 ' incorporates the 
relativistic effects. We define p as the "density" of entities so that p dV gives the 
number of entities within the cross-hatched region in Figure la. (The density p is 
proportional to the fraction of the shell's surface that emits gamma rays during the 
shell's evolution, and therefore gives some idea of how efficiently the shell converts 
its kinetic energy into gamma rays.) The envelope is [5] 

V(T) = if T < T 

T~>a+3 T^a+3 

= V Ta+1 ° if T < T < T max (5) 

T^a+3 T^a+3 

= i> maX Ta+1 ° if T max < T < T 2 (l + P) 2 T 



Time (in units oi T.) 

FIGURE 2. The volume seen by the detector per AT is shown as the solid line, and the 
resulting photon flux at the detector is shown as the dashed line. The volume, T(T), and the 
detector signal, V(T), both increase as T 2 for T < T < T max ; this is identical to the results for a 
non-relativistic shell. However, in the case of a relativistic shell for T > T max , the volume remains 
constant over a large range of T, and the signal shows a power-law decay. Since the detector sees 
a constant volume, and thus a constant number of entities, per AT for T > T max , the decrease in 
the photon flux after T max is entirely determined by relativistic effects. 

where ip is a constant (see Figure 2). 

Equation 5 gives the average signal that would be expected from a collection 
of entities scattered over the surface of a relativistic shell. However, the number 
of entities seen in any volume T(T)dT is a random quantity; since the entities 
are independent and discrete, the randomness follows simple Poisson statistics. 
Therefore, this model predicts time histories with a mean given by Equation 5 and 
Poisson variations about that mean, which would look very similar to the "peaks" 
usually observed in GRB time histories. A low value for p gives a low number 
of entities and correspondingly large Poisson fluctuations, leading to a "spiky" 
time history with many peaks. Conversely, a high efficiency leads to a smooth 
time history. To illustrate these points, we present two simulated time histories in 
Figures 3a and 4a with characteristics selected to match particular BATSE bursts. 
We estimate that the typical fraction of a single relativistic shell's surface that 
emits gamma rays during its lifetime is ~ 10~ 3 [6]. 



CONCLUSIONS 

In order to achieve the observed variability in GRBs using a single-shell model, 
we have found that the gamma-ray emission must occur from small patches on the 
shell (entities). We have derived two significant characteristics of the time histories 
expected from such a model. First, the average number of entities contributing 
to the signal remains constant throughout much of the time history, although the 
overall photon flux decreases due to relativistic effects. Second, the "peaks" in a 
time history can be ascribed to Poisson variations in the actual number of entities 
contributing to the signal at any given time. Taken together, these properties imply 
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FIGURE 3. a. (left) Simulated time history using a high density of entities. Nearly 100% of 
the shell's surface emitted gamma rays in this simulation. Note that the high efficiency gives a 
smooth profile, b. (right) BATSE trigger 1467. 
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FIGURE 4. a. (left) Simulated time history assuming a low density of entities. Only 1% of the 
shell's surface emitted gamma rays in this simulation. Note that this low efficiency gives a spiky 
profile, with many "peaks", b. (right) BATSE trigger 678. 

that the relative variations in a GRB (i.e. the heights of the peaks relative to the 
average envelope of the signal) should remain constant throughout a time history. 
Qualitatively, this result is consistent with visual inspections of GRBs. In general, 
bursts are equally "spiky" during the first half of the time history as the second 
half. 
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